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1. Introduction 

The most plausible candidate for a quintessence explanation of dark energy is a pseudo- 
Nambu-Goldstone boson (PNGB) of a spontaneously and explicitly broken global U{1) 
symmetry Denoting the scale of the explicit breaking by M and the scale of the 

spontaneous breaking by /, the effective mass of the PNGB is rries ~ M"^ / f . The required 
value of this mass for a successful quintessence model is of order the Hubble parameter 
-f^o ~ 10"'^^ eV; this can be obtained naturally if the spontaneous breaking scale / is 
very large, roughly comparable to the Planck scale Mpianck — 10^^ GeV. Because of the 
symmetry, the PNGB has only derivative couplings to matter and radiation, plus couplings 
whose dimensionless strength is suppressed by powers of M// ~ 10~^^. This evades a 
number of strong experimental and observational constraints on weakly-coupled ultralight 
scalars 

In a more general class of PNGB quintessence models, the effective mass of the PNGB 
will vary over time, i.e. it will be a function of the scale factor a obtained by solving the 
coupled cosmological equations of motion. Since / is comparable to the Planck scale, it 
is possible that mefr was much larger at early times, without losing the key property that 
^es/f ^ 1- This raises the possibility that the quintessence PNGB field may also have 
been responsible for primordial inflation, albeit in some modification of the usual scenario. 

A simple avenue towards quintessential inflation is then to assume that mes and the 
PNGB scalar potential V scale (at least roughly) like some power of the Hubble rate H. 
Any model with nics ~ H has the additional virtue of removing the coincidence problem, 
i.e. the fact that the ratio of scales M/f is of order Hq is no longer a coincidence, but 
rather has some dynamical origin. 
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Any model of quintessential inflation must explain why the energy density of the uni- 
verse was dominated by radiation at the time of Big Bang Nucleosynthesis (BBN), even 
though the inflaton dominates the energy density now and dominated it in a primordial 
epoch as well. For a PNGB, the simplest explanation is that the PNGB decays to matter 
and radiation via derivative couplings of the form 

^ijJT^ 9''"d^^^n-^j ' (1-1) 

jWpianck 

where 9 denotes the PNGB field rescaled by / to make it dimensionless, and the Xij are 
moderately small dimensionless couplings. For / ~ Mpianck; this will allow large entropy 
production from PNGB decays. Since the PNGB potential varies over time, the equation 
of state of the PNGB also varies. Thus it is natural to have periods of inflation interspersed 
with periods of radiation dominance. 

The matter coupling ( |l.l| ) was introduced by Cohen and Kaplan in their thermody- 
namic model for baryogenesis Q. In their scenario B or L violating processes occur via 
dimension six four-fermion operators suppressed by a relatively low scale A ~ 10® GeV. 
Combined with a PNGB possessing dimension five couplings like ( |1.1D , they generate a 
baryon asymmetry in thermal equilibrium (at temperatures ^ 10® GeV), and a further 
asymmetry at lower temperatures from PNGB decays. 

In this paper we construct a natural model of PNGB quintessential inflation that also 
implements the Cohen-Kaplan mechanism for baryogenesis. At the same time that we 
provide quintessence and inflation, our model gives a simpler explanation of the baryon 
asymmetry than the original scenario of We are able to assume that the scale for 
the dimension six B, L violating operators is of order 10^^ GeV, or of order 10^^ GeV 
for purely L violating operators. Thus we are slightly above or saturating the current 
experimental bounds from proton decay [^. Our cosmological evolution begins with generic 
initial conditions, unlike the original models of |||]|2|, at an initial temperature Tq ~ 10^^ 
GeV, which is comfortably less than / ~ Mpianck- Our model requires no unnatural tunings 
other than that of the cosmological constant, a tuning that is unavoidable since quintessence 
does not solve the cosmological constant problem. 



2. FRW cosmology driven by a Nambu-Goldstone boson 

Consider a theory with a U{1) global symmetry under which some complex scalar field 
$(x,t) transforms as 

$^e*"$, (2.1) 

where a is a constant. Other fields, including fermions, may also transform nontrivially; 
in particular the symmetry may be a chiral symmetry. We imagine that this symmetry 
is spontaneously broken at some high scale, near Mpianckj determined by the vev / of the 
scalar. Expanding around this vev gives 

$ = (/ + cT(x,i))e^^(^'*) , (2.2) 



- 2 - 



where o" is a real scalar and 9 is another real scalar which has been rescaled by 1// to be 
dimensionless. The effective theory at lower energies has the original symmetry nonlinearly 
realized, with the Nambu-Goldstone boson (NGB) 9 undergoing a shift: 

9{x,t) ^9{x,t) + a . (2.3) 

Obviously the NGB only has derivative couplings in this effective theory. Higher dimension 
operators are suppressed by powers of a high scale / which we are roughly equating to 
-^Planck- To leading order in / the NGB action is just 

d^x^]^fg^^d^9d,9. (2.4) 

Let us further suppose that the energy density of the universe during some epoch is domi- 
nated by the NGB. Then the Priedmann equation is approximately 

= ^9' , (2.5) 
where k'^ = -^pianck/^^- '^^^ equation of motion for 9 (again to leading order in /) is just 

9 + 3H9 = 0, (2.6) 
while the continuity equation is given by 

pe = -6Hpe . (2.7) 

The cosmological solutions for the spatially averaged vacuum expectation value of 9{t) 
depend upon the initial conditions. If we assume that 9 is initially zero, with 9 taking an 
arbitrary initial value, then the solution to (|2.5| )-( |2?7| ) is a static universe. If instead we 
assume that both 9 and 9 have arbitrary initial values, we get an expansion dominated 
(initially) by kination: 



9{t) = -lna{t) , H = ^^, a{t)=t"\ (2.8) 

where we have taken / = ^/^k to simplify notation. Of course in this solution po dilutes like 
1/a^, as appropriate for kination, i.e., an equation of state dominated by kinetic energy. 
Derivative couplings of 9 to ordinary matter will allow the 9 vacuum energy to be converted 
to a thermal radiation bath via decays. Since ordinary matter and radiation dilute like 
1/a? and they will eventually dominate the expansion. 

Let us interpret this solution in terms of the original global symmetry. The equation 
of motion for 9 is just the statement that the global current is covariantly conserved; the 
corresponding conserved global charge of the vacuum Q is proportional to a? 9, which from 
( |2.8| ) is indeed seen to be a constant. Once we include matter couplings the vacuum charge 
Q is no longer constant, but the total global charge including matter contributions is. 
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3. FRW cosmology driven by a pseudo-Nambu-Goldstone boson 



We can modify the discussion in the previous section by introducing a nonvanishing po- 
tential for the 9 field and a noncanonical kinetic function: 

"1 
2' 



-F(cos 9) gf"'d^9du9 - F (cos 9) 



(3.1) 



F and V explicitly break the global U{1) symmetry down to a discrete periodic remnant: 

9^9 + 2-kN , (3.2) 

where is any integer. This kind of explicit breaking would arise if terms proportional to 
powers of <I> + <I>* were present in the original action. Alternatively, there could be Yukawa 
couplings or derivative couplings of $ to fermions: 

x^J<^>^p^iJj, xy^^'d^^i^a.i^j (3.3) 

In this case chiral symmetry breaking will induce, at the loop level, an effective action of 
the form (|3.1| ). The simplest possibility, considered in ||l|, ||, gives 

^1 



d X \f—g 



^fg'"'df,9d,9-M''{l-cos9) 



(3.4) 



9 + 3H9 + -rrrsin 9 = . (3.5) 



where M is some chiral symmetry breaking scale much smaller than /. 

Because of this explicit breaking, 9 is now a pseudo-Nambu-Goldstone boson, with a 
mass of order M"^ / f . It is technically natural for this mass to be small, since it vanishes 
in a symmetry limit. 

Now consider FRW cosmology driven by such a PNGB. The PNGB equation of motion 

is: 

The cosmological solutions again depend upon the initial conditions. One possibility, 
considered in |l], §], is that 9 is initially zero, with 9 taking an arbitrary initial value. 
Then the equation of motion (^]^) is heavily overdamped, and 9 remains approximately 
constant until H(t) decreases to the point where H ~ M'^/f. Thus the PNGB behaves 
approximately like dark energy. The coincidence problem is not solved unless one finds a 
rationale for why M^//, a ratio of two seemingly independent scales, is roughly equal to 
the Hubble rate today. Even then one also needs to explain why the initial value of 9 is 
much less than H. 

Another possibility is that both 9 and 9 have generic initial values during some epoch 
where the PNGB dominates the expansion. Then for M/f <^ 1 the cosmological solution 
is an oscillatory perturbation of the NGB solution ( ^.8]) . The global charge of the vacuum 
Q (X a'^ 9 is no longer conserved; in fact it increases with time like 

This case is more generic than the first one, but does not provide an explanation of dark 
energy. 
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4. Prom generic PNGBs to Slinky 



The most general effective action for the PNGB invariant under the shift (|3.2| ) is 



d^x \/—g 



^ F(cos 9) P{X) - y (cos 9) + Cr, 



(4.1) 



where X = g^^^ d^9dy9, F and V are arbitrary functions, and P is an arbitrary polynomial. 
The matter lagrangian C could contain both derivative couplings to 9 and couplings to 
functions of cos 9. If there are further explicit breakings of the global symmetry, the action 
may also contain a nonperiodic dependence on 9. 

The Slinky model of quintessential inflation introduced in [0-||8[ is a PNGB quintessence 
model of the general class just described, but with some special features. It does not require 
tuned initial conditions; in particular the initial value of 9 is of order H. It is the simplest 
model with a periodic equation of state parameter for quintessence. The quintessence en- 
ergy density dominates the Friedmann equation during an earlier epoch as well as during 
the present epoch, causing primordial inflation as well as present-day acceleration. In order 
for this to happen, both the potential and the kinetic energy must be proportional to the 
square of the Hubble parameter times periodic functions: 

V (xH^ , 9^ (xH'^ . (4.2) 

The first relation implies that these models will have a PNGB mass that decreases with 
time. 

The explicit form of the Slinky is easily derived. We begin with the general action 



(4.1) and make the simplest nontrivial choice P{X) = X. Without loss of generality, we 
will take 9{t) = 1 today. For quintessence the equation of state parameter w^t) should be 
close to —1 today; for simplicity we will take it to be exactly —1. Since 

l + w{t) = F{cos9)9^ . (4.3) 

This means that F{cos 9)9^ ^ now. The simplest nonvanishing function that does this 
is 

F = ^(l-cos0), (4.4) 
where 6 is a dimensionless constant given by 

3 Afpianck 




/ 



(4.5) 



The explicit breaking represented by the cosine term in ( ^.41 ) does not have any tun- 
able small parameter associated with it. Thus one might worry that this breaking is not 
small. However since this breaking occurs in the kinetic function, it generates global charge 
nonconserving processes that are suppressed by powers of momenta divided by k. Thus for 
dynamics well below the Planck scale this breaking is indeed small. 
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The PNGB potential V{6) must vanish when we turn off the explicit breaking. Of 
course this requirement is a tuning, reflecting the fact that quintessence models do not 
explain why the ^-independent cosmological constant vanishes. In this limit we are kination 
dominated, with w{t) = —cos 9 = 1. Thus the simplest form for the potential is 

V o^{l + cos9)H'^{e) . (4.6) 

We need to solve three equations of motion, beginning with the Friedmann equation: 

= ^P9 . (4.7) 

The second equation is the scalar EOM: 

o = e + 3He + ^ye^ + ^ , (4.8) 

where a prime denotes a derivative with respect to 6. The third equation is the continuity 
equation: 

pe = -3(1 + w)Hpe , (4.9) 

where w{9) is the equation of state parameter for quintessence. 
Plugging in ([4.4|) and (|4.6|) we get a solution with 



H = Ho el;^^-^''^^) (4.10) 



where Hq is the Hubble parameter today. 

Now we have enough information to specify F and V completely: 



F{e) = _(l-cos0). 
Vie) = i/?o(l + cos0)exp 



1(0- sin 0) 





(4.11) 



where po = 3k'^HQ/2 is the analog of in the standard PNGB quintessence model 
discussed earlier. Note that the scalar potential has an additional explicit breaking that 
does not respect the periodicity ( |3.2| ). This breaking is small during any epoch such that 
H/k < 1. 

Slinky is a PNGB quintessence model characterized by a periodic equation of state 
w = —cos6. The effective mass of the PNGB is tied to the expansion rate, and therefore 
decreases over time. These two features combined allow the same PNGB to serve the 
role both of the inflaton in the early universe and the quintessence field that drives the 
accelerated expansion in the current epoch. 

The number of inflationary epochs that have occurred is controlled by the dimensionless 
parameter b; larger values of b mean more inflationary periods. Most values of b are ruled 
out by the requirement that the universe had at most a tiny inflationary component during 
BBN [|. Note from (|4|) that b < 0.49 implies that / > Mpianck- The Shnky models 
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Figure 1: The cosmological history of the simplest SUnky model. The x-axis is the logarithm base 
10 of the scale factor. Shown are energy densities of radiation (green), dark matter (orange), baryons 
(blue) and dark energy (red), as a fraction of the critical density. 



considered in have values in the range 0.09 < b < 0.4. This does not necessarily mean 
that these models require trans-Planckian vevs; for example, if we replace the single PNGB 
by N identical ones, then b is effectively rescaled as 6 ^ h/\fN . Thus small values for b do 
not indicate a scenario that is technically out of control. 

The nonstandard cosmological history of this model is depicted in Figure |l|. The energy 
densities of radiation, dark matter, baryons and dark energy, computed as a fraction of 
the critical density, are plotted as function of the logarithm of the scale factor. We have 
chosen to start the cosmological evolution at a scale factor of lO"^'^, but this choice is not 
essential. We have assumed that the Friedmann equation is initially dominated by the 
PNGB. For simplicity we set the initial radiation density to zero, so all radiation arises 
from PNGB decays, as explained in the next section. 

The radiation temperature as a function of the scale factor is shown in Figure ^. During 
the initial inflation the temperature is roughly constant at around 10^^ GeV, increasing 
mildly towards the end of this first inflationary epoch. This behavior is in contrast to 
standard primordial inflation, where the temperature first decreases very rapidly, then 
increases very rapidly (reheating). The new behavior is due to the time- varying equation 
of state. 



5. Couplings to matter 



The allowed couplings of 9 to matter can be divided into two classes: 
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Figure 2: The cosmological history of the simplest SUnky model. The x-axis is the logarithm base 
10 of the scale factor. Shown are energy densities of radiation (green), dark matter (orange), baryons 
(blue) and dark energy (red), as a fraction of the critical density. 



1. Couplings of functions of cos 6 to matter. These couplings explicitly break the global 
U{1) symmetry. It is technically natural to take the dimensionless coupling constants 
for all such couplings to be small. 

2. Derivative couplings of to matter. Some of these couplings respect the full global 
C/(l) symmetry, so there is no symmetry argument for tuning their dimensionless 
coupling constants to be small; indeed we will assume that they are of order one. 
However almost all such couplings are higher dimension operators suppressed by 
powers of momenta over powers of Mpianck- 

Thus for Slinky, as well as for more generic PNGB quintessence models, it is a good 
approximation to only include couplings to matter of the form 

Xij g^"'d^,e^a^^j , (5.1) 



where Xij is dimensionless. A caveat is that the Slinky potential (4^) has an additional 
explicit breaking of the global U{1) symmetry, that we are assuming has no analog in the 
matter couplings. 



Couplings of the form ( p.lD allow the 6 field to decay into ordinary matter. This 



process can be modeled, albeit roughly |10|, as an additional friction term in the 9 equation 
of motion: 

o = e + 3He + Te + l^e^ + ^ . (5.2) 

2 r r 
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Here T is the decay width, which for decays into pairs of fermions can be genericahy written 

r = kmruefiit) , (5.3) 

where km ~ Xij are dimensionless couphngs which we will take to be approximately con- 
stant with magnitudes in the range .1 to .01; rries oc H is the time varying mass of the 
PNGB, obtained by expanding the potential ( |4.6| ). 

This process converts quintessence vacuum energy into matter and radiation. Thus 
Slinky produces the following nonstandard cosmological history: 

• During a primordial epoch, w is close to -1 and inflation occurs. 

• Eventually w increases towards +1, and the radiation produced from 6 decays, which 
is now diluting less rapidly than the quintessence energy, comes to dominate the 
expansion. 

• The process repeats. In the original Slinky model [^, which we will use from now 
on, the second radiation dominated epoch overlaps with the time of Big Bang Nucle- 
osynthesis. 

• Well before BBN time, nies becomes so small that the PNGB decays to matter effec- 
tively turn themselves off, due to kinematic suppression. PNGB decays to photons 
continue, but these are loop-suppressed and so do not give large entropy production. 

• Right now we are entering the third inflationary phase. 

From the arguments given above, the only couplings in this model whose values affect 
the cosmological evolutions are b and the Xij. The value of b is coarsely adjusted such that 
we ensure that BBN time is radiation dominated and that w is close to —1 today. The 
dominant Xij couplings are adjusted such that the radiation and matter fractions Qr/^A 
and Qdm/^a come out to their measured values today. For simplicity we assume that 
all of the dark matter is produced thermally, as e.g. in standard WIMP scenarios. The 
remaining matter couplings (for the moment) remain free. 

We now see more clearly why the the temperature history represented in Figures || and 
|2| is so nonstandard. The coupling between matter and the PNGB field force them track 
each other, giving significant entropy production even after the first period of inflation . 
This entropy production turns off well before BBN time due to the kinematic suppression 
of PNGB decays. 



6. Baryogenesis 

If baryon number B and lepton number L are exactly conserved, then B or L violating 
decays of the PNGB simply do not occur. This can be seen from the coupling (5T), which 
with integration by parts vanishes for any conserved current = '07^ V'- 



-9- 



However we certainly do not expect the global symmetries B and L to be respected by 
unification scale physics. We therefore introduce all possible B and L violating dimension 
six operators: 

suppressed by some superheavy scales Aij^i- Suppressing flavor labels, this implies that 
the B or L violating decay width of the PNGB can be written 

r = ^ . (6.2) 

Thus PNGB evolution and decay during the first inflationary epoch will produce a net 
B and L asymmetry, including a, B — L asymmetry that will survive the B + L washout 
by sphalerons during the electroweak phase transition. 

At sufficiently high temperatures the B oi L violating processes will be in thermal 
equilibrium. The equilibrium condition is [Q] 

T:^i = ^{l-cose)HpNGB , (6.3) 

where Hpi\fGB denotes the Hubble rate that would result from the PNGB alone, given in 
equation (4.10). Figure (0) shows an estimate of the value of the scale factor at which B 



or L violating processes decouple, assuming A = 10^^ GeV. The corresponding decoupling 
temperature is about 8 x 10^^ GeV. 

Once we are out of equilibrium the baryon asymmetry continues to grow via PNGB 
decays. The net baryon number density as a function of time can be written as |Q] 

n^ = x!I^}^esme. (6.4) 
It is useful to express the baryon number density as a function of 9, i.e. 

«B = A ^ -p- sm^ . (6.5) 

To calculate the baryon asymmetry produced, we have to include the expansion of the 
universe. This can be easily done by replacing all the volume factors in the equations above 
by the corresponding comoving volumes, proportional to a? , 

- ia^ns) = \a^j^-^sm9. (6.6) 

To obtain the net baryon number density we have to numerically integrate this expression 
from the decoupling temperature to 9f, the temperature where baryon production from 
PNGB decays is cut off kinematically by the decrease in rueff- For decays to baryons 9f 
corresponds to a scale factor of approximately 10~^^, when meg ~ H{9) ~ 1 GeV. 
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Figure 3: The ratio of the relative rates of B or L violating processes versus the cosmological 
evolution of the PNGB field, as a function of the logarithm base 10 of the scale factor. We have 
taken A — 10^^ GeV. The ratio is shown in red; B or L violation goes out of equilibrium below the 
blue line. 

From the baryon density we want to extract the baryon to photon ratio, rj. As usual 
the photon number density is given by 

n, = ^ (6.7) 

Contrary to the standard scenarios [0, 12], where the baryon (more precisely, B — L) to 
photon ratio does not change after baryogenesis, in our model there is significant entropy 
production due to the coupling of the PNGB field to radiation. The net effect of such 
a production will be to dilute any baryon to photon ratio produced before BBN, where 
entropy production stops. Therefore r/ can still be calculated at Of, but we have to include 
the extra dilution factor 7, given by 

7 = ■ (6.8) 

Here S = T^, with the effective number of relativistic degrees of freedom. Alter- 
natively, we can calculate 77 directly at BBN time. 

The only input parameters left to adjust in our simple model are the scale A and the 
dimensionless coupling A. We consider two well-motivated scenarios. In the first scenario, 
A ~ 10^^ GeV, just about saturating the generic lower bounds from the nonobservation of 
proton decay |^. We estimate the baryon to photon ratio n] in this case to be 

?7 ~ 4A X 10"^° . (6.9) 
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This gives the observed baryon asymmetry for A ~ 1. 

The second scenario has A ~ 10^^ GeV, as might be appropriate in a leptogenesis 
model. Here we estimate 



V 



5A X 10"^ , (6.10) 



which gives the observed baryon asymmetry for A ~ 0.01. 

The cosmological history of r] is unhke that of previously discussed models of quintessen- 



tial baryogenesis [^]-|2C]. At early times the PNGB effective mass is large, and the baryon 
number violating processes are least suppressed. A baryon or lepton asymmetry is pro- 
duced in thermal equilibrium, but this contribution is negligible compared to the asymme- 
try produced later by PNGB decays. This is basically because the PNGB decay process is 
enhanced by the square of the Planckian scalar vev, as seen in (^^). The net baryon 
production turns itself off kinematically as rrieff decreases. The baryon to photon ratio 
drops dramatically during the subsequent second inflation era, where there is further large 
entropy production. The fact that the baryon to photon number density observed today is 
very small is thus due entirely to the existence of a second inflation era; B or L violation 
was only mildly suppressed at the time that most of the net baryon or lepton excess was 
created. 

7. Conclusions 

We have shown that a single noncanonical PNGB could be responsible for both primordial 
inflation and the present day accelerated expansion, while simultaneously generating the 
observed baryon excess. The baryon asymmetry is generated after the first inflationary 
period, from PNGB decays via dimension six operators that violate B and/or L. The ratio 
of baryon to photon number densities is greatly diluted later on, via entropy production 
from B and L conserving PNGB decays. The baryon energy density is completely negligible 
expect during two eras: the present day and an era around the time of the electroweak 
phase transition. 

While economical, technically natural, and consistent with current data, our model has 
some theoretical shortcomings. It explains dark energy but does not solve the cosmological 
constant problem. There is also no first principles explanation for why the PNGB potential 
scales like the square of the Hubble rate. Resolving these shortcomings presumably involves 
an ultraviolet completion of the PNGB effective theory into a more fundamental framework. 

One of the distinguishing features of this model, is that it predicts substantial entropy 
production in the era between the electroweak phase transition and BBN. If dark matter is 
predominately composed of thermally produced WIMPs, this prediction can be tested by 
combining collider data with signals from direct and indirect WIMP searches Q, pl[|-p5|. 

Because of the entropy production after the second inflation era (but before BBN), 
any pre-existing baryon asymmetry is much diluted. This requires that the operators 
responsible for B violation be not too much suppressed, forcing us close to the current 
experimental bounds for proton decay. Thus another prediction is that proton decay will 



be observed in one of the future proposed experiments |26| 
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Future dark energy probes will pin down the equation of state of dark energy with 
much greater precision. In our simple model we artificially set w = 1 exactly at redshift 
z = 0; more relevant is that w varies by about 2% as redshift is varied between and 2. 
This variation is comparable to the one-sigma projected combined errors after the Stage 
IV dark energy probes |30|. 



Models of the type discussed here predict running of the spectral indices of the Cosmic 
Microwave Background (CMB) |^T|. These effects may be large enough to extract from 
observations of the PLANCK satellite [^]. The smoking gun of a predictive model of 
quintessential inflation is that these CMB effects, a result of inflation, are directly related 
to the detailed equation of state of dark energy. 
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